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ABSTRACT 

We present Hubble Space Telescope optical coronagraphic polarization imaging observations of the 
dusty debris disk HD 61005. The scattered light intensity image and polarization structure reveal a 
highly inclined disk with a clear asymmetric, swept back component, suggestive of significant inter- 
action with the ambient interstellar medium. The combination of our new data with the published 
1.1 Ilia discovery image shows that the grains are blue scattering with no strong color gradient as 
a function of radius, implying predominantly sub- micron sized grains. We investigate possible ex- 
planations that could account for the observed swept back, asymmetric morphology. Previous work 
has suggested that HD 61005 may be interacting with a cold, unusually dense interstellar cloud. 
However, limits on the intervening interstellar gas column density from an optical spectrum of HD 
61005 in the Na I D lines render this possibility unlikely. Instead, HD 61005 may be embedded in 
a more typical warm, low-density cloud that introduces secular perturbations to dust grain orbits. 
This mechanism can significantly distort the ensemble disk structure within a typical cloud crossing 
time. For a counterintuitive relative fiow direction — parallel to the disk midplane — we find that the 
structures generated by these distortions can very roughly approximate the HD 61005 morphology. 
Future observational studies constraining the direction of the relative interstellar medium fiow will 
thus provide an important constraint for future modeling. Independent of the interpretation for HD 
61005, we expect that interstellar gas drag likely plays a role in producing asymmetries observed in 
other debris disk systems, such as HD 15115 and 5 Velorum. 

Subject headings: circumstellar matter — planetary systems: formation — planetary systems: proto- 
planetary disks — stars: individual (HD 61005) 
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1. INTRODUCTION 

Nearly two dozen dusty debris disks surrounding 
nearby stars have now been spatially resolved at one or 
more wavelengths. Many of these systems show clear 
similarities. For example, the radial architecture of sev- 
eral debris disks can be understood in terms of a uni- 
fied model of steady-state dust production via collisions 
in a parent planetesimal belt (e.g., Strubbe & Chiang 
2006). However, while t he observed structure of many 
systems is ring-like (Kalas et aO2006l : IWvattll2008f ) . most 
disks show substructure such as clumps, warps, offsets, 
and brightness asymmetries not explained in traditional 
steady-state collisional grinding models. 

These unexpected features have triggered a great deal 
of recent theoretical work. The effects of massive plan- 
etesimal collisions, sandblasting by interstellar grains, 
close stellar fiybys, dust avalanches, and secular and 
resonant perturbations by exoplanets have all been in- 
voked to explain the observations (e.g., Moro-Martin et 
al. 2007, and references therein). However, as many of 
these theories produce similar structures, further obser- 
vational constraints are needed to better understand the 
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key forces affecting disk structure and the circumstances 
in which they apply. 

At a heliocentric distance of 34.5 pc (iPerrvman et al.l 
I1997D . the debris disk surrounding HD 61005 (SpT: G8 
V; Gray et al. 2006), is a promising target for advanc- 
ing our understanding in this area. The significant in- 
frared excess for this source (Lir/L* = 2 x 10~^) was 
recently discovered as part of the Spitzer FEPS sur- 
vey (Carpenter et al. 2009), indicating 60 K blackbody- 
emitting grains ^16 AU from the star. Follow-up Hubble 
Space Telescope (HST) coronagraphic imaging observa- 
tions with the Near Infrared Camera and Multi-Object 
Spectrometer (NICMOS; HST/GO program 10527; D. 
Hines, PI) resolved the source (Hines et al. 2007, here- 
after H07), revealing an unprecedented swept, asymmet- 
ric morphology, suggestive of significant interaction with 
the interstellar medium (ISM). H07 suggested that this 
system could be a highly inclined debris disk, undergo- 
ing ram pressure stripping by the ambient ISM. However, 
this interpretation requires an unusually high interstel- 
lar density for the low-density Local Bubble in which HD 
61005 resides. Furthermore, the single wavelength inten- 
sity image was insufficient to provide strong constraints 
on the dominant size of the scattering grains and the 
overall scattering geometry. 

To further quantify the physical properties of grains 
seen in scattered light and the overall geometry of the 
system, we obtained optical coronagraphic polarimetry 
imaging observations of HD 61005 with the Advanced 
Camera f or Surveys (ACS) on board HST. As demon- 
strated bv lGraham et al.] (|2007l ) for the case of AU Mic, 
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polarization observations in scattered light are invalu- 
able for breaking degeneracies between grain scattering 
properties and their spatial distribution. Furthermore, 
the ACS data represent a factor of two improvement in 
angular resolution compared to the 1.1 jiia discovery ob- 
servations. In addition to these new imaging data, we 
also obtained a high resolution optical spectrum to char- 
acterize ambient interstellar gas surrounding this system. 
In §2, we describe the steps taken in observing and re- 
ducing these data. In §3, we discuss the results of these 
observations, their consequences for the system scatter- 
ing geometry, and the additional constraints they pro- 
vide when combined with the 1.1 /im NICMOS image. 
In §4, we explore whether interactions with ambient in- 
terstellar gas can plausibly explain the observed swept, 
asymmetric morphology in this system. We discuss the 
implications for these potential explanations in §5 and 
summarize our findings in §6. 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained optical coronagraphic observations of HD 
61005 using the ACS high resolution camera (HRC) I'.'S 
diameter occulting spot on 2006 December 19 (HST/GO 
program 10847; D. Hines, PI). In each of two contiguous 
orbits, we imaged HD 61005 with the F606W filter in 
combination with the POLOV, POL60V and POL120V 
polarizer filters (two 340 second exposures per filter com- 
bination). Aside from the telescope position angle, which 
is rotated 23.032 degrees between orbits, the observa- 
tional procedures for the two HD 61005 orbits were iden- 
tical. Before the HD 61005 orbits, we observed the point- 
spread function (PSF) reference star HD 82943 (SpT: F9 
V; Gray et al. 2006; B = 7.16, V = 6.56) using an observ- 
ing sequence identical to the two orbits allocated to HD 
61005 (B = 8.93, V = 8.22). We also observed a second 
PSF star, HD 117176 (SpT: G5 V; Gray et al. 2001; B 
= 5.69, V = 5.00), in an identical manner following the 
HD 61005 orbits. 

For each filter combination, we combined the two 340 
second, pipeline processed (bias subtracted, fiatfielded) 
frames by excluding the maximum value at each pixel 
position, thereby minimizing the impact of cosmic ray 
events. After dividing by the cumulative integration time 
of each frame, we performed sky subtraction by taking 
the median value in a 10x20 pixel box in the lower left 
corner of the chip, which is the position farthest from 
the bright target star. We registered the images by se- 
lecting a fiducial HD 61005 image (the POLOV image 
in the first orbit) and subtracting all other HD 61005 
frames using small offsets (0.02 pixels) to minimize the 
residuals in regions dominated by light from the stellar 
PSF. The offsets that minimize residual differences be- 
tween frames were then applied to the individual POLOV, 
POL60V, and POL120V images to align them to a com- 
mon reference frame relative to the star. We carried out 
an identical registration procedure for the two PSF ref- 
erence stars. 

We then subtracted the HD 61005 PSF in each of the 
three POLOV, POL60V, and POL120V frames by the 
corresponding frames from each of the two PSF refer- 
ence stars. Prior to subtraction, we scaled each refer- 
ence star to match the expected brightness of HD 61005, 
using photometry obtained from the direct images. The 
HD 117176 observations, made immediately following the 



HD 61005 orbits yielded a better subtraction than the 
HD 82943 observations, made five weeks prior to the HD 
61005 orbits. We therefore used the subtraction obtained 
with HD 117176 for all subsequent analysis. 

Following PSF subtraction, we corrected the resul- 
tant images for geometric distortion yielding 25 mas x 
25 mas pixels. We then constructed Stokes parameter 
images corrected for instrumental polarization following 
lPavlovskvl([2006h . For the ACS/HRC/F606W/POLV in- 
strumental configuration and high fractional linear polar- 
ization {p = (Q^ + U'^Y^'^/I > 0.2), the residual system- 
atic error is 10% of the computed polarization fraction. 
For less strongly polarized sources {p < 0.2), the sys- 
tematic error in the degree of linear polarization is ap- 
proximately constant at ap = 0.01. In both cases, the 
systematic uncertainty in position angle is 3°. 

We next calculated polarization vectors from the 
derived Stokes images. As the polarization frac- 
tion is intrinsically positive and biased upwards by 
noise, we employed the spatial binning algorithm of 
iCappellari fc CopinI (|2003l ) to bin the Stokes /, Q, and U 
images to approximately constant signal-to-noise prior to 
this calculation. Within ~ 1.4^' of the star, the computed 
polarization vectors become significantly disordered in 
magnitude and direction, as a result of systematic PSF 
subtraction errors. We therefore only consider polariza- 
tion vectors outside this radius in our analysis. The sur- 
face brightness at 0.9 — 1.4'' is similarly compromised and 
should be treated with caution. 

Finally, we converted from instrumental to physical 
brightness units using the synthetic photometry pack- 
age, Synphot. As input, Synphot requires the instrument 
configuration (camera, coronagraph, wideband and po- 
larizing filters) and the source spectrum across the band. 
Since the latter is unknown, we performed the calcula- 
tions three times assuming: (1) a fiat spectrum, (2) a 
Teff = 5500 K Kurucz synthetic spectrum approximating 
the G-dwarf stellar spectrum, and (3) a spectral slope 
across each band that is the same as that calculated be- 
tween the NICMOS and ACS bands using method (2). 
All methods yielded conversion factors within 1% of each 
other, suggesting the assumed source spectrum factors 
negligibly into the total color uncertainties. 

In addition to the ACS observations, we analyze two 
additional data sets: (1) The NICMOS FllOW image; 
a full description of the NICMOS data acquisition and 
reduction is given in H07. (2) High resolution (R ^ 
60,000) echelle spectra for HD 61005 and two compari- 
son stars of similar spectral type (HD 33822: Teff = 5850 
K, HD 13836: Teff = 5580 K; Masana et al. 2006). The 
spectra were obtained on the Keck I telescope with the 
HIRES spectrometer on 30 Dec 2004, 16 Jan 2006, and 
09 Nov 2008. The wavelength range was 3700-6200 A, 
though our analysis concerns only the Na I D lines at 
5889.951 A and 5895.924 A. We used a standard proce- 
dure to perform flat-fielding, sky subtraction, order ex- 
tr act ion , and wavelength calibration of the raw echelle 
images (jButler et al. I[l996l : l\^t et al.llT99l . 

3. RESULTS 

3.1. ACS Scattered Light and Polarization 

3.1.1. Two-dimensional morphology and polarization 
structure 
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Fig. 1. — Logarithmically-scaled image of HD 61005 in Stokes / 
at > 0.1 mJy arcsec"^ with polarization vectors and a 1.8^^ coro- 
nagraphic mask overplotted. The color bar units were calculated 
using the synthetic photometry package, Synphot. The plotted 
polarization vectors were computed from Stokes images binned to 
approximately constant signal-to-noise. 

Figures [T] and [2] display the F606W total intensity im- 
age of HD 61005 on a logarithmic and linear scale, re- 
spectively. The figures show two distinct morphological 
components. The first component, denoted by NEl and 
SWl in Figure [21 resembles a near edge-on disk. The 
putative midplane for this component is observed to ex- 
tend out to ~ 3'' from the star, where the signal-to-noise 
per pixel falls below unity. The second morphological 
component is detected below the nominal disk midplane. 
This component exhibits an asymmetric, "swept back" 
morphology, suggestive of significant interaction with the 
interstellar medium. This unusual structure was first 
noted by HOT in their NICMOS FllOW discovery image. 
Both the disk-like and swept back components addition- 
ally exhibit a striking asymmetry between the northeast 
and southwest sides of the source. At a given projected 
radius, the northeast side of the source is approximately 
twice as bright as the southwest side. This brightness 
asymmetry is also seen in the NICMOS FllOW image 
(e.g., HOT, Figure 4). 

Figure [3] displays the polarized fiux of image of HD 
61005, while Figure [1] overplots polarization vectors on 
the Stokes / image, spatially binned according to the 
procedure outlined in §2. The polarization vectors in 
Figure [T] show that along the plane of the disk, the frac- 
tional polarization increases with radial distance from 
the star from ~ 10% to ~ 35%. The orientation of the 
electric field within ~ 2.1'' also appears perpendicular 
to the nominal disk midplane. In the swept back area 
of the source, the polarization vectors rotate to become 
approximately perpendicular to the outer edge of this 
component. The fractional polarization in this region is 
similar to that seen in the outer part of the disk compo- 
nent. 

3.1.2. Disk component parameters 

The ACS polarization results add further evidence to 
the suggestion by HOT that HD 61005 is a near edge- 
on debris disk. In particular, the HD 61005 polariza- 
tion structure is very similar to that seen in the spec- 







70 AU 


NEl 






NE2 










\^SW1 


N 
t 




\sW2 


E <J 






4 


2 

6RA (arcsec) 


-2 -4 



Fig. 2. — Same as Figure [T] only with a linear display scale and 
without polarization vectors. Overplotted are labels for the surface 
brightness components displayed in Figure [5] 
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Fig. 3. — Logarithmically-scaled polarized flux {\/ Q"^ + U"^) im- 
age of HD 61005. 

trally blue edge-on debris disk around the M dwarf, AU 
Mic, observed using the same instrumental configura- 
tion CGraham et al.1 120071 ). The midplanes of both HD 
61005 and AU Mic exhibit high fractional polarization 
(Pmax ^ 0.4) which increases with projected radius. Sim- 
ilarly, both disks exhibit an electric field orientation per- 
pendicular to the disk midplane at all projected radii. 
These features are expected for small-particle {x < 1) 
scattering in an edge-on disk (Kruegel 2003). 

For HD 61005, these effects are quantified in Figure 
m which shows the binned polarization vector position 
angles and magnitudes from Figure [TJ Linear fits to the 
one-dimensional polarization position angles (PA) versus 
projected distance within 2.1^' imply fiat slopes, as ex- 
pected for a highly-inclined disk geometry (—5.1 ± 2.5 
and 1.0 ± 3.3 for the northeast and southwest sides, re- 
spectively). Averaging all polarization position angles 
within 2.1^', the implied disk position angle is 71.7 ± 0.7 
deg, where the listed uncertainty is purely statistical and 
does not include the additional 3 deg calibration uncer- 
tainty (§2). The steady increase in polarization fraction 
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Fig. 4. — One dimensional version of the binned polarization 
vectors displayed in Figure [T] In the top panel, the data show a 
clear increase in fractional polarization with distance from the star. 
Linear fits to the data are overplotted in grey with a solid line for 
the southwest side and a dot-dashed line for the northeast side. In 
the bottom panel, the rotation of the electric field orientation in 
the swept back outer disk is seen in the deviation from the nominal 
midplane polarization position angle outside ~ 2.1" . Linear fits for 
data in and outside this radial distance are overplotted. 



with projected radius is also seen in Figure [H Linear fits 
to the binned polarization fraction over the full extent of 
the source imply slopes of 0.09 ± 0.01 and 0.11 ± 0.01 for 
the northeast and southwest disk sides, respectively. 

As an independent check on our disk interpretation, 
we fit elliptical isophotes to the total intensity image 
from Figures [T] and [21 assuming a circularly symmet- 
ric disk viewed in projection. We independently fit eight 
isophotes outside 1.4'' for / > 1.3 mJy arcsec"^. The 
implied disk position angle from these fits is 70.7 ± 0.5 
deg, in good agreement with the position angle inferred 
from the electric field orientation within 2.1" . The im- 
plied inclination to the line of sight from these fits is 
i = 80.3 ± 0.6 deg. 

3.1.3. Swept component parameters 

The disk structure described above contrasts with pre- 
viously imaged presumed interstellar dust phenomena, 
such as the infrared bow structure surrounding the A 
star 6 Velorum (Caspar et al. 2008) and the filamentary 
cirrus surround ing some Vega-lik e stars with significant 
infrared excess (jKalas et al.ll2002[ ). Nevertheless, Figures 
[T] and [2] clearly reveal a second asymmetric component 
of the source not typical of nearby debris disks and sug- 
gestive of interaction with the interstellar medium. 

The polarization signature of this component is evident 
in the systematic rotation of the polarization vectors out- 
side ~ 2.1''. Linear fits to the one-dimensional polariza- 
tion position angles versus projected distance (Figured]) 
give non-zero slopes of 8.2 dz 1.0 for the northeast lobe 
—5.8 ±1.9 for the southwest lobe. Estimates for the po- 
sition angles of the outer edge of this component (cuts 
NE2 and SW2 in Figure [2]) are obtained from the out- 
ermost polarization vector position angles. The implied 
position angle for component NE2 is 82.8 ± 2.1 deg; the 
result for SW2 is 65.1 =b 3.8 deg. 



3.1.4. One- dimensional surface brightness profiles 

One-dimensional surface brightness profiles through 
the disk midplane (components NEl and SWl in Fig- 
ure [21 PA = 70.7 deg) and along the swept component 
outer edge (components NE2 and SW2 in Figure [21 PA = 
82.8 deg, 65.1 deg) are shown in Figure [5l The midplane 
surface brightness follows a broken power law. The fitted 
power law indices between 0.9" and 2.1" are — 1.1 ± 0.1 
for both midplane disk lobes. Between 2.1" and 3.5", 
power law fits give indices of —7.2 ± 0.3 for component 
NEl and —6.0 ± 0.2 for component SWl. Breaks in the 
scattered-light surface brightness of debris disks are of- 
ten taken as the location of the parent body population 
for the scattering grains (e.g., Strubbe & Chiang 2006). 
However, as the midplane break approximately coincides 
with the position at which the outer edge of the swept 
component intersects the disk, this standard interpreta- 
tion may not hold for this case. 

The surface brightness profile for the outer edge of the 
swept component is well represented by a single power- 
law (Figure [5]). The fitted power law index is — 4.1 ± 0.1 
for both NE2 and SW2. The listed errors for the power 
law fits in this section are formal fitting errors and should 
be treated as lower limits to the true uncertainties. 

3.2. ACS+NICMOS 
3.2.1. Disk Scattered Light Colors 

We computed the color of the disk by rebinning the 
ACS Stokes / image to the same pixel resolution as 
the NICMOS image (0.0759"). We then convolved the 
binned ACS image with a coronagraphically unocculted 
field star (the approximate NICMOS PSF) and per- 
formed the corresponding operation on the NICMOS im- 
age. Finally, we divided each image by the stellar flux 
density of HD 61005 at the appropriate band-center ef- 
fective wavelength before dividing the NICMOS image 
by the ACS image. 

Figure [6] illustrates that the HD 61005 debris disk ap- 
pears predominantly blue with no significant systematic 
color gradient. The mean intensity ratio inferred from 
Figure [6] is 0.32 ± 0.10 (corresponding to a color index 
of [F606W]-[F110W] = -1.2 ± 0.3), where the error is 
dominated by uncertainties in the PSF subtraction. Al- 
though comparing images of different resolution can re- 
sult in systematic color errors (e.g., Colimowski et al. 
2006), our convolution steps appear to have a small ef- 
fect; the mean color neglecting convolution is within la 
of that found including convolution. 

The blue color inferred from Figure [6] (ratio < 1) is 
rare, as the handful of debris disks with color measure- 
ments to date show mainly red colors (Meyer et al. 2007, 
and references therein). There are several notable excep- 
tions, however. The HD 32297 and HD 15115 debris 
disks, for example, have been suggested to s how blue op- 
tical to near-infra red scattered light colors (|Kalas|[2QQ5l : 
iKalas et al.ll2007[ ). though t he resu l t for H D 32297 is cur- 
rently under debate (Debes et al." 2009). Interestingly, 
both disks have morphological features consistent with 
ISM interaction. HD 32297 shows a bowed disk struc- 
ture, similar to HD 61005, though on a much larger 
scale (- 1000 AU; Kalas et al. 2005). HD 15115 is 
highly asymmetric, perhaps as the result of ISM erosion 
(■Debes et al....2009l 



HD 61005 Optical Imaging 



5 




Fig. 5. — Surface brightness profiles through the disk midplane (left) and the outer edge of the swept back component (right). The cuts 
are one pixel wide, and the components listed in the legends are labeled in Figure [2] The directions for the cuts are taken from the position 
angles derived in §3.1.2 and 3.1.3. The power law fits described in the text are overplotted in grey. The data show a pronounced brightness 
asymmetry between the northeast and southwest disk lobes. 
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Fig. 6. — ACS and NICMOS images processed to compute the grain color. The left panel displays the ACS image binned to the NICMOS 
pixel resolution and convolved with the NICMOS off-spot PSF. The middle panel displays the NICMOS image convolved with the ACS 
PSF. The right panel shows a masked ratio image of the left and middle images divided by the stellar fiux density ratio; values less than 
unity represent grains that preferentially scatter blue light, whereas values greater than unity represent grains that preferentially scatter 
red light. In the ratio image, we have applied a mask to all pixels with values less than 2.5 times the background level in either original 
image. The ratio image indicates that the disk appears predominantly blue with no appreciable color gradient. 



The M dwarf debris disk, AU Mic, which has a 
similar polarization structure to HD 61005 (§3.1.2), 
also shows blue optical to near-infrared colors, with a 
color gradient towards bluer colors at larger radial dis- 
tances, indicating change s in the ^rain size distr ibution 
(IStrubbe Chiand l2006l : iFitzgerald etHI [2007h . The 
global [F606W]-[F110W] color of the HD 61005 disk is 
comparable to the [F606W] — J AU Mic disk color at pro- 
jected radii within ~ 40 AU, the approximate location 
of th e parent body ring in this system (Fitzgerald et al. 
[20071) . 

Like AU Mic, the blue color of the HD 61005 disk is 
likely due to the disk grain size distribution. Scattered 
light images afford a relatively narrow window on the 
grain size population because the brightness at a given 
location in the disk is largely determined by the prod- 
uct of the grain size distribution and the grain scatter- 
ing cross section. In the presence of a steep size spec- 



trum characteristic of a collisional cascade, this product 
is strongly peaked near x = ~ 1. For example, for 

the D ohnanyi spectrum with dn/da ex a~^/^ (iDoh nanvil 
Il969l ). this peak occurs at x :^ 2 — 6, depending on the op- 
tical properties of the grain material. In Figure [71 a Mie 
calculation shows the peak for water ice and astronomi- 
cal silicates at 0.6 /im and 1.1 /im as a function of grain 
size. The plots illustrate that the HD 61005 NICMOS 
images trace grains with radii of order 0.2 — 2 /im, while 
the ACS images trace systematically smaller grains with 
radii of order 0.1 — 1 /im. Thus the observation that HD 
61005 is globally brighter at optical wavelengths than 
near-infrared wavelengths suggests the disk contains a 
larger number of grains at increasingly small sizes, con- 
sistent with our expectation of a steep size spectrum. 

The inference that the F606W ACS images trace pre- 
dominantly sub-micron sized grains is also consistent 
with the imaging polarimetry results (§3.1). The large 
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Fig. 7. — Product of the scattering cross section and a Dohnanyi 
size spectrum as a function of logarithmic grains radius for spher- 
ical particles. At a given wavelength, this product largely deter- 
mines the surface brightness at a particular location in the disk. 
Assuming a steep size spectrum (e.g., Dohnanyi), a large peak is 
observed near x ~ 1. Thus a given scattered light image traces 
a relatively narrow window of the grain size population. The HD 
61005 disk is globally brighter at optical wavelengths than near- 
infrared wavelengths, likely reflecting the larger number of grains 
at increasingly small sizes. The above results are monochromatic: 
^F606W = 0-6 l-J^rn and Apiiow = 1-1 M"^- 

polarization fraction and electric field orientation per- 
pendicular to the edge-on midplane are in qualitative 
agreement with the expected signature of scattering by 
small spherical particles with x < 1 (Kruegel 2003). For 
larger spherical grains, the electric field orientation can 
rotate by 90° at certain scattering angles, resulting in an 
orientation parallel to the edge-on midplane. Further- 
more, any line of sight comprising emission from a range 
of scattering angles will tend to show weak linear polar- 
ization (Grah am et al.l [20071 ). Neither of these features 
is consistent with the HD 61005 polarization results. 

By integrating the curves in Figure [3 and comparing 
the results to the measured color, we can constrain the 
size distribution. The results of this procedure are shown 
in Figure [HI which shows the implied color for ice and 
silicate grains for grain size distributions of the form: 
dn/da (X ^ where amin < a < 1 mm. A range of min- 
imum grain sizes are considered, as radiation pressure 
could potentially remove a fraction of small grains from 
the system (Appendix A) . The data place an upper limit 
on the minimum grain size: amin ^ 0.3 /im. Figure[8]also 
suggests the global size distribution is steeper than the 
Dohnanyi spectrum: a ~ 4.5 — 5.5. This slope is consis- 
tent with results from collisional equilibrium modeling of 
other debris disk systems (e.g., Strubbe & Chiang 2006). 
However, the applicability of these models to the unique 
HD 61005 system is presently uncertain. 

3.2.2. Astrometric test for low-mass companions 

The ACS and NICMOS images contain four field stars 
in common to both data sets, roughly offset from HD 
61005 in RA and Dec by [-2A" , 9 A" ], [-8.6'', 0.0'' ], 
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Fig. 8. — Implied disk colors for ice and silicate grains for grain 
size distributions of the form: dn/da oc a~^, where amin < a < 
1 mm. The data place an upper limit on the minimum grain size 
(«min ^0-3 /im) and suggest a global size distribution steeper than 
the canonical Dohnanyi size spectrum (a = 3.5). 

[-12.8", 3.5"], and [-14.0", 2.4"]. As the ACS and first- 
epoch NICMOS images were obtained 1.1 years apart, 
we can measure the field star proper motions to investi- 
gate whether any of these stars are Ukely companions to 
HD 61005. The annular proper motion of HD 61005 is 
fia = —56.09 mas yr~^ in right ascension and /j^s = 74.53 
mas yr~^ in declination (Ferryman et al. 1997). Thus co- 
moving companions are expected to show a 102.61 mas 
displacement between the two epochs of observation. 

To derive the NICMOS stellar positions, we first cal- 
culated stellar centroid positions on the image frames 
uncorrected for geometric distortion using the apphot 
center task, as described in Co x et al.l (1997). We then 
applied a distortion correction to the raw positions using 
the correction coefficients appropriate for Cycle 15 data, 
given in the NICMOS data handbook (version 7.0). The 
reported rms astrometric uncertainty from applying pre- 
viously derived distorti on corre ctions to early commis- 
sioning data is 13.6 mas (I Cox et al. 1997). Data obtained 
more recently yield similar results (e.g., Schneider et al. 
2006). 

The geometric distortion for the ACS data re- 
quires a low spatial frequency correction for the op- 
tical telescope assembly (OTA) and ACS optics, and 
two high frequency corrections for the given wide 
band filter and polarizing filter (lAnderson fc Kin^ 
12004 Kozhurina-Platais & Biretta 120041 ). Follow- 
ing [Kozhurina-Platais fc Biretta (^20041 ). we first used 
the effective PSF libr ary and fitting technique of 
"Anderson & King (2000) to derive raw positions for the 
field stars from the flat-fielded images {^-f It. fits) ob- 
served through the POLOV filter. We nex t applied the so- 
lution derived bv Anderson fc Ki ng* ("2004) to obtain stel- 
lar positions corrected for the low frequency OTA distor- 
tion and the high frequency F606W filter distortion. Fi- 
nally, we applied a further correction for the distortion in- 
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tra duced by the PQLOV filter, us i n^ the solution derived 
by iKozhurina-Platais fc Birettal (|2004f ). The reported 
rms precision derived from applying; this method to com- 
missi oning data is 1 mas (jKozhurina-Platais fc Birettal 
[200l . 

To test for possible companionship, we performed rela- 
tive astrometry by adopting one star's position as a fixed 
reference point and calculating the relative proper mo- 
tion for the remaining three stars. We repeated this pro- 
cedure three times, using each field star as the reference. 
No systematic motion for any of the field stars is ob- 
served; all relative proper motions are less than twice the 
approximate expected rms positional accuracy (cr = 13.6 
mas). We note that although the distortion solutions 
employed only strictly apply to noncoronagraphic, di- 
rect imaging data, the small measured astrometric offsets 
suggest the additional field distortions imposed by the 
coronagraphic optics are negligible in this case. Given 
that the annual proper motion of HD 61005 is signifi- 
cantly greater than the measured astrometric motions of 
all field stars, we conclude that unless all four sources are 
companions, the four field stars are background objects. 

4. INTERPRETATION 

The results presented in HOT and the previous section 
strongly suggest that HD 61005 is a near edge-on de- 
bris disk, undergoing significant erosion by the ambient 
interstellar medium. In this section, we explore whether 
interaction between the disk and local interstellar gas can 
plausibly explain the observed swept, asymmetric mor- 
phology. 

4.1. Interaction with a Cold, Dense Cloud 
4.1.1. Ram pressure stripping of bound grains 

In their discovery paper, HOT suggested that interac- 
tion with a cold (T ~ 20 K), dense (n ~ 100 cm"^) cloud 
could potentially explain the HD 61005 morphology. In 
such a cloud, ram pressure on disk grains from inter- 
stellar gas could unbind grains from the system, analo- 
gous to the process t hat strips g:as from cluster galaxies 
(jGunn fc Got^lT972l : Ivan Gorkoml[200l . 

For ram pressure stripping to operate, the drag force 
on a grain must be comparable to or greater than the 
gravitational force binding the grain to the star. For a 
grain of radius agrain and density Pgrain orbiting at dis- 
tance r from a star of mass Mgtar^, the interstellar cloud 
density n and relative cloud-disk velocity v must obey: 



/ n \ / V 
V200 cm-3y V30 kms-i 



> 



a 



agrain 



0.95 Mo/ VO-lMm 
70 AU 



2 cm 



(1) 



This required density is characteristic of cold, dense 
gas, the existence of which is constrained in §4.1.3. In 
such a high-density cloud, HOT argued that Bondi-Hoyle- 
Lyttleton (BHL) accretion could also potentially play a 

^ Throughout this paper, we adopt Mstar = O.95M0, 
based on the pre- m ain-se quence evolutionary tracks of 
ID'Antona Mazzitellj (pWzl ) and (Baraffe et al. (1998) and 
the FEPS age estimate reported by HOT (E. Mamejek, private 
communication) . 



role, leading to an accumulation of interstellar grains that 
contribute non-negligibly to the observed infrared excess 
emission and scattered light morphology. 

4.1.2. Ram pressure deflection of unbound grains 

Recently, iDebes et al.l ()200 9') suggested that ram pres- 
sure deflection of unbound grains could plausibly shape 
several previously resolved debris disks, including HD 
61005. Their model does not specify the origin of 
the unbound population. However, such a substantial 
population of unbound grains are unlikely to be pro- 
duced in a steady state situation, as collisional equi- 
librium models predict that the scattered light surface 
brightness due to bound grains dominat es over th at 
from grains unbound by r adiation pressure (|Krivov et all 
120061 : iStrubbe fc Chian'3 [20061 Furthermore, for the 
case of HD 61005, radiation pressure from the low lu- 
minosity star may be insufficient to unbind grains of any 
size (Appendix A). 

Let us nonetheless suppose that a substantial un- 
bound grain population exists. The ambient ISM density 
needed to deflect such grains by the distances implied 
by our observations is similar to that required to strip 
bound grains ( §4.1.ip . Both scenarios require densities 
characteristic of cold, dense clouds. For /3 = 1 grains 
launched from parent bodies on circular orbits at 70 AU, 
the required interstellar cloud density and relative veloc- 
ity obey^: 



( ^ ^ ( ^ Y > ( 

V100cm-3y V20kms-i/ ^ \0.95 



95 Mo 

Pgrain 

2gcm~^ 
y/x_ 
0.2 



'-f'gram 

0.1 /im 
70 AU ^ 

110 AU 



X 



(2) 



Here x and respectively, are the relative distances trav- 
eled by the grain parallel and perpendicular to the disk 
midplane at a given time after the grain is born. We 
adopt X ~ 110 AU and y ~ 20 AU, the approximate 
positions of the outermost binned polarization vectors in 
Figure [H 

4.1.3. Limits on cold, dense interstellar gas 



The explanations of §4.1.11 and §4.1.21 are hampered 
by two factors. First, in both circumstances, escaping 
grains leave the system on orbital timescales (10^ — 10^ 
yr). These timescales are shorter than the timescale over 
which the disk is expected to be collisionally replenished 
with sub-micron sized grains (> 10^ yr; see Appendix 
B). Thus both scenarios require that we are observing 
HD 61005 during a short-lived period in its history. 

Second, cold neutral medium (CNM) clouds within the 
Local Bubble are extre mely rare, occupyi ng; a volumetric 
ffihng factor of <C 1% (Welsh et al.'1994). To date, only 
one cloud with properties typical of CNM clouds (n ~ 50 

^ Note that equations [T] and [2] assume the cross section grains 
present to interstellar gas equals the geometric cross section; cal- 
culations performed using ballistic cluster-cluster aggregates and 
ballistic particle-cluster aggregates suggest this assum ption is good 
to within an order of magnitude (,Minato et al.ll20Q6h . 
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Fig. 9. — Optical spectrum of HD 61005 and two comparison late- type G star spectra in the Na I D lines. The Na I lines in HD 61005 
are very similar to the comparison star lines. No narrow absorption component expected for an intervening dense cold interstellar cloud 
is observed, placing limits on the potential effectiveness of ram pressure stripping of bound grains and deflection of unbound grains in the 
HD 61005 system. 



T - 20 K, V = 4.9 X 1.4 x 0.07 pc^) is known 
within the Local Bubble {d < 45 pc, towards 3C 225; 
Heiles et al. 2003; Meyer et al. 2006; Meyer 2007). Fur- 
thermore, in contrast to optical stellar spectra towards 
this known cloud, our optical spectrum of HD 61005 does 
not show evidence of an analogous CNM cloud. This 
finding is illustrated in Figure [H which shows spectra 
of HD 61005 and two comparison late-type G stars (HD 
33822, HD 13836) in the Na I D lines (D2: 5889.95 A, 
Di: 5895.92 A). HD 61005 is a relatively young star with 
a detected rotational velocity (Vsinz = 9 km s~^; Holm- 
berg et al. 2008). Thus for direct comparison purposes, 
we have convolved the spectra of both HD 33822 and 
HD 13836 with a Gaussian of this width. In Figure [9l 
we have also effectively removed the stellar Doppler shift 
of the comparison star spectra by cross-correlating each 
spectrum with the HD 61005 spectrum and shifting it by 
the appropriate amount. 

Figure [§] shows that the HD 61005 spectrum is very 
similar to the late-G comparison star spectra in the 
Na I lines. No narrow absorption component, char- 
acteristic of the known CNM Local Bubble cloud to- 
wards 3C 225, is observed (see Figure 1 in Meyer 2007). 
The corresponding 2a upper limit to the Na I gas col- 
umn is log(A^Nai/cm~^) < 10.6, based on the contin- 
uum signal-to- noise ratio (S/N ~ 100), the instrumen- 
tal resolution (R ~ 60,000), a linewidth characteristic 
of the nearby cold cloud towards 3C 225 {b = 0.54 km 
s~^), and assuming the linear regime of the curve of 
growth. This upper limit translates into a total hydro- 
gen column of log(A^Hi+H2/cni~^) < 19.0, employing the 
sodium-hydrogen conversion relation derived by Ferlet 
et al. (1985). The scatter in this relation is such that a 
third of measured hydrogen columns deviate from their 
predi cted columns by factors of a few ( Wakker fc MathisI 
l2000f ). Taking this into account, our predicted hydrogen 
upper limit is still well below the columns expected for 
CNM clouds (10^0 - 10^1 cm-2; McKee & Ostriker 1977) 
and also below the mean column density of Na I detected 
towards 3C 225: (log(iVNai/cm-2)) = 11.7 (.Mever et al., 
l2006h . 



There are a few examples of tiny (< 1000 AU) low col- 
umn density (log(A^Hi+H2/cm~^) < 19.0) cold clouds in 
the literature (e.g., Heiles 1997, Stanimirovic & Heiles 
2005). These clouds constitute < 2 - 4% of the total 
neutral hydrogen column along the lines of sight in which 
they are detected. The precise frequency of these clouds 
is not well constrained by large-scale surveys, as their 
low columns are similar to survey detection thresholds 
(- 10^^ cm-2; Heiles & Troland 2005). However, an up- 
per limit to their frequency can be obtained by noting 
that > 5 such clouds along a typical line of sight would 
lead to their systematic detection in large-scale surveys 
(C. Heiles, private communication). Thus within the Lo- 
cal Bubble (d < 100 pc), these tiny, dense clouds occupy 
< 0.025% of the distance along a typical line of sight. At 
present, only two dozen debris disks have been spatially 
resolved. Therefore it seems unlikely that we have al- 
ready imaged a debris disk residing in one of these clouds. 
Interaction with a cold, dense cloud is thus an unsatis- 
factory explanation for the swept morphology observed 
in HD 61005. 

4.2. Interaction with a Warm, Low-Density Cloud 

Given that disk interaction with a cold, dense cloud 
appears unlikely, we explore in this section whether in- 
teraction with a warm, low-density cloud can potentially 
explain the observed morphology. 

4.2.1. Ram pressure stripping of disk gas and entrainment 
of disk grains 

Warm (T ~ 7000 K) interstellar clouds dominate the 
mass of the local interstellar medium (e.g., Frisch 2004) 
and occupy a local volum e filling factor of ~5.5%— 19% 
(|Redfield fc Linskvl |2008[ ). Our optical spectrum does 
not constrain the presence of such a cloud towards HD 
61005, as typical columns towards these clouds are Nm ~ 
10^^ cm-2 (Redfield 2006). 

The low densities of warm local clouds are insufficient 
to supply a ram pressure force comparable to the gravita- 
tional force and thereby directly unbind grains from the 
system. For typical parameters of warm local interstel- 
lar clouds (nm = 0.2 cm~^, Vrei = 25 km s~^; Redfield 
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2006) and 0.1 /im disk grains, the ram pressure stripping 
radius is ~ 10^ AU (Equation 1), much larger than the 
radius of the observed bow structure (~ 2" = 69 AU). 

However, in principle, ram pressure stripping could still 
play a role if the HD 61005 disk contains gas that is 
undergoing ram pressure stripping by the ISM. In this 
scenario — in contradistinction to the direct ram pres- 
sure stripping scenario outlined in §4.1.11 — disk grains are 
swept away by the interstellar flow only because they are 
entrained in outflowing disk gas. Gas-gas interactions 
correctly explain the HI morphologies of galaxies travel- 
ling through an intracluster medium (van Gorkom 2004); 
the truncated, swept-back HI disks of galaxies strongly 
resemble HD 61005 (e.g., see Figures 1.7 and 1.8 of van 
Gorkom 2004). We show below, however, that this inter- 
pretation for HD 61005 is incorrect because the require- 
ment on the density of disk gas is incompatible with the 
requirement that grains be entrained. 

To unbind disk gas, the ISM ram pressure must exceed 
the gravitational force per unit disk area: 

nfiv^ > GMstar^/r' (3) 

where 1^ ^2x 10-^^ g is the mean molecular weight of 
the ISM and a is the surface mass density of disk gas. 
For parameters appropriate to a warm cloud (iRedfieldl 
[2006), all disk gas having a surface density 




is swept away by the ISM. The circumstellar gas content 
of HD 61005 is unknown. 

For unbound disk gas to entrain disk grains, the mo- 
mentum stopping time of a grain in gas cannot be 
much longer than the outflow timescale, 1/1^, over which 
marginally unbound gas departs the system, where Q 
is the Kepler orbital frequency. From Weidenschilling 
(1977), the momentum stopping time of a grain in rar- 
efied gas is given by the Epstein (free molecular drag) 
law as 

1 flgrainPgrain 500 / ^^grain \ f Pgrain \ 




which is too long compared to the outflow time. Thus 
we discount the possibility that the observed disk mor- 
phology arises from ram pressure stripping of disk gas. 

4.2.2. Secular perturbations to grain orbits induced by ram 

pressure 

The direct ram pressure stripping scenario described in 
§4.1.11 requires that the ISM ram pressure force on disk 
grains be comparable to the stellar gravitational force. 
As shown in §4.2. H this condition is not met in low- 
density warm clouds, for which Fram/^grav ^ 10~^ (at ~ 
70 AU; see Equation 1). However, even in the case when 
the force exerted by interstellar gas is much less than 
the gravitational force, neutral gas can introduce secular 



perturbations to bound grain orbits that could signifi- 
cantly change the morphology of the disk over timescales 
of ~ 10^ — 10^ yr, assuming sub-micron sized grains dom- 
inate the scattered-light distribution. This perturbation 
timescale is less than the crossing time of local warm 
clouds : tcross ^ (Lcioud/5pc)/(vrei/25kms~-^) ~ 10^ yr 
(jRedfie ld 2006). Thus this mechanism can plausibly ex- 
plain the disturbed HD 61005 morphology. Interestingly, 
this scenario has been proposed as the primary removal 
mechanism for dust in our own solar system at 20—100 
AU (Scherer 2000.), though at present, little empirical 
evidence is availa ble to test this theory. 

As described in^Sc hererl (|200Ql ). the underlying physical 
process responsible for neutral gas drag is similar to that 
responsible for solar wind drag. In both cases, momen- 
tum transfer from incident protons or gas particles to 
the grain surface results in secular perturbations to the 
grain's initial orbit. However, the monodirectional char- 
acter of the interstellar gas drag force leads to changes in 
particle orbits that are very different from those induced 
by the solar wind. While the solar wind and Poynting- 
Robertson drag both act to reduce grain eccentricities 
and semi-major axes, interstellar gas drag tends to in- 
crease them, eventually unbinding the grains from the 
system. 

In the abse nce of o t her p erturbing forces, the ana- 
lytic work of [Scherer (2000) shows that the gas drag 
force acts to rotate a given particle's orbital plane into a 
plane coplanar to the flow vector, and its star-pericenter 
(Runge-Lenz) vector into a direction perpendicular to 
the flow vector. These effects are illustrated in Figures 
[TOl and [m which show the orbital evolution of two dust 
grains with different initial orbital elements. The grain 
in Figure [TOl starts on a low-eccentricity (e = 0.3) orbit in 
the X — y plane with its Runge-Lenz vector anti-aligned 
with the incoming interstellar flow and its angular mo- 
mentum vector aligned with the z— axis (out of the page). 
This orientation causes the grain to be decelerated on one 
leg (e.g., point pi) and accelerated on the other (e.g., 
point P2), causing the Runge-Lenz vector to rotate. The 
grain's periastron at each successive time is denoted by 
a cross to highlight this effect. The rotation continues 
until the Runge-Lenz vector becomes perpendicular to 
the incoming flow. Thus, counterintuitively, neutral gas 
drag leads to a build up of grains perpendicular to the 
relative flow direction. 

The tendency of neutral gas drag to rotate a given 
grain's orbital plane into a direction coplanar with the 
incoming flow is illustrated in Figure [TTJ which shows the 
evolution of a low-eccentricity (e = 0.1) grain initially in- 
clined 80° with respect to the ISM flow. The grain's ini- 
tial Runge-Lenz vector is chosen to be at its equilibrium 
orientation (perpendicular to the flow) to isolate the ef- 
fect of the orbital plane rotation. The ISM flow exerts 
a net torque on the orbiting grain, causing the angular 
momentum vector to rotate into a direction perpendicu- 
lar to the incoming flow. This effect is easiest to discern 
from the two-dimensional projection of the y — z plane 
in Figure [m 

The above discussion suggests that rigorously model- 
ing the neutral gas drag effect requires knowledge of the 
initial grain orbital elements and the interstellar gas flow. 
However, for the case of HD 61005, neither of these pre- 
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Fig. 10. — Orbital evolution of a test particle in the x — y plane 
under the effect of neutral gas drag. Five time periods separated 
by a fixed number of orbits are shown. The grain begins with its 
star-pericenter (Runge-Lenz) vector anti-aligned with the incoming 
interstellar flow, directed along the positive y— axis, and its angu- 
lar momentum vector aligned with the z— axis (out of the page). 
The ISM flow decelerates the grain on one leg (e.g., point pi) and 
accelerates it on the other (e.g., point ^2), causing the Runge-Lenz 
vector to rotate. The grain's periastron at each successive time is 
shown by a cross to illustrate this effect. 




Fig. 11. — Orbital evolution due to neutral gas drag of a test 
particle initially inclined by 80° with respect to the incoming flow. 
Five time periods separated by a flxed number of orbits are shown. 
The grain's initial Runge-Lenz vector is chosen to be at its equilib- 
rium orientation. The ISM flow exerts a net torque on the orbiting 
grain, causing the angular momentum vector to rotate into a di- 
rection perpendicular to the ISM flow. 



requisites is known. HD 61005 has a well-determined 
space motion; the Hipparcos-measured proper motion 
corresponds to a plane-of-sky velocity of = 9.2 ± 0.3 
km s~^, vs = 12.2 ± 0.4 km s~^ at the distance of 
HD 61005. The radial velocity is Vr = 22.3 ± 0.2 km 
s~^ ([Nor dstrom et al. 2004). However, the velocity of 
the putative cloud responsible for the swept morphol- 
ogy is unknown. Velocities of local warm clouds can be 
comparable to the observed space motion of HD 61005 
(|Frisch et al.l[2002h . Thus the red arrow in Figure 3 of 
H07 denoting the direction of the star's tangential mo- 
tion is not a reliable indicator of the cloud-star relative 
velocity. 

Unfortunately, HD 61005 is difficult to assign to any 
known interstellar clouds, owing to its heliocentric dis- 
tance (34.5 pc) and galactic coordinates (/ = 246.4°, b = 
-5.6°). Redfi eld fc Linskvl (2008) recently used radial 
velocity measurements of 157 lines of sight to identify 15 
warm clouds in the local interstellar medium. However, 
the identified clouds are thought to reside largely within 
15 pc of the Sun, whereas the distance to HD 61005 is 
34.5 pc. Furthermore, the line of sight to HD 61005 is 
not assigned to any of these clouds. The star's galactic 
coordinates could plausibly associate it with either the 
G cloud or the Blue cloud (e.g., see Figure 19 in Redfield 
& Linsky 2008), though HD 61005 is more likely asso- 
ciated with a more distant, currently unidentified cloud, 
rendering the cloud-star relative velocity for this system 
highly uncertain. 

The initial orbital elements prior to the presumed cloud 
interaction are similarly unknown. The radial architec- 
tures of debris disks in scattered light show significant di- 
versity and currently appear largely indepen dent of other 
observables, such as age or spectral t ype ()Kalas et al.l 
^06). Furthermore, as discussed in §3.1.4[ the one- 
dimensional surface brightness profiles of HD 61005 do 
not allow us to place strong constraints on the location 
of the parent body population, and thus the grain eccen- 
tricities and inclinations. We do note, however, that the 
sub-micron particle sizes implied by the polarization and 
color measurements do not a priori contradict the bound 
orbits requirement of the gas drag perturbation theory, 
as radiation pressure from the star is likely insufficient 
to unbind grains of any size (Appendix A). 

Given the large uncertainties in the disk and ambi- 
ent ISM properties of the HD 61005 system, we are un- 
able to empirically test whether gas drag is responsible 
for the asymmetric, swept morphology. Therefore to ex- 
plore whether this mechanism can plausibly explain the 
observed str ucture, we ado pt the numerical techniques 
described in Scherer (2000) to construct a modest grid 
of models with reasonable assumptions for the disk grains 
and ISM. Details of the model construction and restric- 
tions are described in Appendix B. In summary, we sub- 
ject a ring of 0.1 /im, low eccentricity (e = 0.2) grains 
with semimajor axes of 60 AU and random inclinations 
to a uniform density cloud typical of nearby warm in- 
terstellar clouds (nm = 0.2 cm~^) traveling at a typical 
cloud-star relative velocity of Vj-ei = 25 km s~^ (Redfieldl 
l2006h . We test a range of relative fiow vectors and pro- 
duce scattered light images from the resultant grain dis- 
tributions after the system has achieved a steady state. 

The resulting grid of models produced for a range of 
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relative flow directions and disk inclinations is presented 
in Figures [12] — [m Each frame is ^" across, and the color 
scale is logarithmic. The adopted coordinate system is 
described in Appendix B. A comparison of Figure [T] to 
Figures [12] — [14] suggests that none of the IS gas drag 
models are a striking match to the data. On the other 
hand, a gross swept morphology, somewhat similar to HD 
61005, is present in a few of the model panels. These 
best approximations correspond to disk inclinations of 
i = ±80°, consistent with the data (§3.1.2), and relative 
flow directions largely coplanar with the disk midplane, 
approximately in the plane of the sky (Figure [121 panels 
0°,180°). This relative flow direction is counter- 
intuitive, going against the flow direction suggested by 
preliminary inspection of the scattered-light images, per- 
haps providing a further strike the against IS gas drag 
interpretation. On the other hand, if these preliminary 
models are a first approximation of the data, then this 
distinctive relative velocity may be verified by future ob- 
servations (§5.1). 

For ease of viewing, we show in Figure [15] the ACS 
and NICMOS data together with a neutral gas drag 
model that roughly approximates the observed morphol- 
ogy, corresponding to the top row, third panel of Figure 
[121 A summary of the morphological shortcomings of 
this strongest IS gas drag model can be summarized as 
follows: 

1. The swept back components NE2 and SW2 are 
much more extended relative to the main disk in 
the data than in the model. 

2. The edges of the swept back structure are more 
pronounced in the data than in the model. 

3. The sharp radial spurs observed beyond the inflec- 
tion point along components NEl and SWl are not 
present in the model. 

4. The model figure does not show the significant 
brightness asymmetry between the northeast and 
southwest disk lobes that is seen in the data. 

These differences could indicate that the physics incor- 
porated into the current IS gas drag models is overly 
simplistic (§5.1), or that interstellar gas drag is not re- 
sponsible for the observed morphology (e.g., §5.3). 

As discussed in §3.1, the ACS polarization results can 
be understood qualitatively through geometric consider- 
ations alone. Thus the most promising IS gas drag mod- 
els naturally reproduce the observed polarization struc- 
ture. 

5. DISCUSSION 

5.1. Interstellar Gas and the HD 61005 Morphology 

The previous section explored whether disk/gas inter- 
action can plausibly explain the unusual HD 61005 mor- 
phology. Of the four scenarios considered, three are im- 
plausible, given the limits on the ambient interstellar gas 
density imposed by our optical spectrum. The fourth 
scenario, secular perturbations from low density gas, is 
an attractive alternative, as this mechanism can signif- 
icantly distort grain orbits well within a cloud crossing 
time. Furthermore, the densities required by this sce- 
nario are typical of local interstellar clouds, which occupy 
up to - 20% of the local ISM. 



Nevertheless, our preliminary modeling of this effect 
(Figures [12] — [14]) can only produce disk morphologies 
in very rough agreement with the observations, suggest- 
ing that either additional physics needs to be incorpo- 
rated into the current models, or that an altogether dis- 
tinct physical mechanism is at work. Indeed, the current 
models are simplistic, and their applicability is limited 
by several key assumptions: 



1. Astrosphere Sizes: As discussed bv iSchererl ( 2000[ ), 
the toy models presented in Figure [12] — [H] re- 
quire that the disturbed grains be inside the as- 
trospheric termination shock, such that the inter- 
stellar gas density and velocity can be approxi- 
mated as constant. For the case of HD 61005, the 
termination shock distance is unknown. Further- 
more, as HD 61005 is farther away than any star 
for wh ich a direct a strospheric detection has been 
made (|Woodl [2004f ), the astrosphere size may be 
difficult to constrain observationally. In general, 
termination shock distances vary greatly, depend- 
ing on the ambient ISM and stellar wind conditions 
(e.g., densities, temperatures, velocities, stellar ac- 
tivity). For the ca se of the Sun only, the hydro- 
dynamic models of iMiiller et al.l ([2006) show that 
the termination shock distance could easily vary 
between ~ 10 AU and 500 AU. Observational as- 
trosphere measurements of solar-type stars are con- 
sistent with these predictions^ (Mann et al. 2006, 
and references therein). 

2. Initial Conditions: As noted in §4.2.2 and Ap- 
pendix B, our adopted initial orbital elements for 
the HD 61005 disk grains prior to the interstel- 
lar cloud interaction are highly uncertain, given a 
lack of information for the grain properties and un- 
derlying planetesimal population that collisionally 
replenishes the observed dust disk. Future long- 
wavelength observations sensitive to larger grains 
may be able to place tighter constraints on the dis- 
tribution of sub-micron grains prior to the inter- 
stellar cloud interaction, as the distribution of large 
grains would likely reflect that of the parent bod- 
ies for the sub-micron population. Furthermore, 
larger grains would not be significantly affected by 
interstellar gas drag on the same timescale as the 
sub-micron size grains traced in these observations. 
The nu n ierica l techniques employed here and in 
IScherei] (|2000h can be easily revised to accommo- 
date an arbitrary initial disk architecture, provided 
the orbits are not highly eccentric, such that aver- 
aging over one orbit and applying Gauss' method 
is invalid. 

3. Internal Disk Collisions: In Appendix B, we esti- 
mate to order-of-magnitude that the collision time 
for submicron grains at ~70 AU is ~5000 yr. 
This collision time is somewhat longer than the 
timescale over which our model relaxes to a steady 
state — essentially the time for gas drag to un- 
bind a grain — given in Appendix B as ~3000 yr. 

^ On a broader scale, no astrosphere detections have been made 
for stars earher than G-type. As a result, the typical effect of 
astrospheres on IS gas drag in debris disks surrounding A-type 
stars is difficult to reliably assess. 
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Fig. 12. — Steady state model images for hypothetical debris disk systems undergoing neutral gas perturbations. In all models, the 
relative flow is coplanar with the disk midplane. The adopted cylindrical coordinate system is shown with respect to the face-on disks in 
the left column. The vector r points in the direction of the relative ISM flow; the azimuthal orientation of r is deflned by 0. The disk 
inclinations are indicated at top. Each box is ^" x ^" (assuming a distance to the system of 34.5 pc); the color scale is logarithmic. The 
models show that brightness asymmetries, bow structures, and swept morphologies can all be produced by disk encounters with warm 
interstellar clouds, which occupy a sizable fraction of the local ISM. 



That the times are comparable supports the as- 
sumption of our models that each grain removed by 
gas drag is collisionally replenished. At the same 
time, the comparison of timescales underscores a 
shortcoming of our model — that removal of grains 
by collisions is ignored. In reality, submicron grains 
should be removed from the system not only by gas 
drag, but also by collisions, in roughly equal pro- 
portions. We defer to future work a comprehensive 
study that includes removal by collisions via a col- 
lisional cascade. 



4. Planeta r y Con figurations: Finally, as illustrated by 
iSchererl (j2QQQl) (e.g., his Figure 2), the incorpora- 
tion of planetary orbits can appreciably change the 
perturbed orbital elements from the case in which 
only IS gas drag is considered. This caveat is es- 
pecially important for massive grains, or grains in 
close proximity to planetary orbits. As a result, 
the models presented here should be treated with 
some caution if applied to typical planetary sys- 
tem scales (< 50 AU; e.g., Kenyon & Bromley 
2004). This cautionary point may be particularly 
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Fig. 13. — Same as Figure [T2] only with components of the relative flow both coplanar with the disk midplane and perpendicular to it. 
The radial and perpendicular flow components are equal in magnitude. 



important for the case of HD 61005, as the origin 
of the brightness asymmetry between the north- 
east and southwest disk lobes (§3.1.1) is unknown. 
The agreement between the northeast and south- 
west deflected component position angles (§3.1.3 
and Figures [1] — [2]) suggests this asymmetry may 
originate from a physical mechanism entirely dis- 
tinct from ISM interaction. If the asymmetry is 
due to a massive perturber, the disk morphologies 
produced in Figures [12] — [HI are likely to be in- 
applicable. Resolved long wavelength observations 
sensitive to massive grains are needed to further 
explore this possibility. 

In addition to the above uncertainties, a remaining 
ambiguity important for future IS gas drag modeling is 



the velocity of the putative cloud responsible for the HD 
61005 morphology. HOT noted that the star's tangential 
space motion is perpendicular to the disk midplane, in 
agreement with the relative flow vector suggested by ini- 
tial inspection of the observed images. However, while 
assigning the relative flow direction to the star's tangen- 
tial velocity is appealing, the three scenarios explored in 
§4 which assume this flow direction were found to be un- 
tenable. Furthermore, velocities of local warm clouds can 
be comparable to the observed space motion of HD 61005 
(§4.2.2). As a result, the star's tangential motion is not 
a reliable indicator of the cloud-star relative velocity. 

Future spectroscopic observations may be able to de- 
tect the cloud directly (e.g., HST/GO Program 11674; 
H. Maness, PI), thereby providing key constraints on the 
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Fig. 14. — Same as Figures [12] and [13] only with the relative flow vector oriented purely perpendicular to the disk midplane. The 
dust distribution in these models is largely axisymmetric, as described in the text. The brightness asymmetries evident at intermediate 
inclinations {i = 45°, 80°, —80°, —45°) are the result of scattering asymmetry; for positive inclinations, the lower half of the disk in the 
image is closer to the observer than the upper half. 




dx (orcsec) 

Fig. 15. — Comparison between the ACS Stokes / image (top), 
the NICMOS 1.1 /im image (middle), and a promising neutral gas 
model (bottom), taken from Figure [12] (top row, third column). 
This comparison shows that the simple model presented here is 
only a very rough representation of the data. Still, the gross swept, 
asymmetric morphology is clearly present in the model image. 

ambient ISM density and velocity. Such observations will 
greatly inform future modeling, as the preliminary inter- 
stellar gas drag models presented here suggest a counter- 
intuitive relative motion parallel to the disk midplane, 
rather than perpendicular to it. Though HD 61005 is dif- 
ficult to assign to any known interstellar clouds (§4.2.2), 
its galactic coordinates could plausibly associate it with 
either the G cloud or the Blue cloud. The space veloc- 
ities of both of these clouds suggest a relative motion 
dominated by the radial velocity component. Therefore 
if HD 61005 is associated with either of these clouds, the 
relative motion is inconsistent with all models posited in 
§4. 

5.2. General Applicability of Interstellar Gas Drag 

The normalcy of the interstellar densities, velocities, 
and cloud sizes required by the secular perturbation 



model in §4.2.2 suggests that IS gas drag can be impor- 
tant beyond HD 61005 in shaping debris disk morpholo- 
gies. Taking the simple models of §4.2.2 and Appendix 
B at face value, several of the general morphological fea- 
tures produced in Figures [T2l — [T4l are consistent with ob- 
served disk structures. For example, the e xtrem e bright- 
ness asymmetry in HD 15115 (jKalas et al.ll2007[ ) may po- 
tentially result from interstellar gas drag, though a range 
of alternative explanations could explain this system as 
well (e.g., see list in §1). The bow structures seen in some 
of the face-on models in Figure [12] are also reminiscent 
of the mid-infrared morphology observed around the A 
star, 6 Velorum, which was rec ently modeled as a p urely 
interstellar dust phenomenon (jGaspar et al.ll20Q8[ ). Fi- 
nally, the middle panels in Figure [13] show that warps, 
similar to that seen in /3 Pic (e.g., Mouillet et al. 1997), 
can in principle be produced for a relatively wide range 
of fiow directions. 

However, while IS gas drag can in principle produce 
commonly observed disk features, the rate at which gas 
drag is expected to affect the observations remains un- 
clear. Beyond the uncertainties in the model physics de- 
scribed in §5.1, the characteristics of warm, low density 
clouds are currently uncertain, as detailed knowledge of 
them is limited to clou ds residing predom inantly within 
15 pc of the Sun (Red field Linskvll2008h . As a result, 
our understanding of typical cloud sizes, shapes, and to- 
tal volumetric filling factor remains rudimentary. A key 
finding in this area, however, is that a significant frac- 
tion of nearby warm clouds appear to exhibit filamentary 
morphologies, which would limit the average interaction 
time between a given disk and cloud, likely reducing the 
rate at which IS gas perturbations produce an observable 
effect. This concern is particularly important for the case 
of disks surrounding early-type stars, as grains traced in 
scattered light tend to be larger in disks surrounding A- 
type stars than in their later type counterparts, owing 
to the larger radiation pressure blowout size. As such, 
the scattered-light morphologies for A-star disks require 
correspondingly longer cloud-disk interaction times to be 
noticeably affected. The timescale for a given grain to be- 
come unbound under IS gas drag incr eases approxi mately 
as the square root of the grain size (|Scherei]l2000h . 

5.3. Interstellar Grains and the HD 61005 Morphology 

Finally, we note that all models posited in §4 con- 
sider only the role of interstellar gas, ignoring the poten- 
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tial ef fects of interstellar grains. lArtvmowicz fc ClampinI 
(|1997l ) investigated IS sandblasting of debris disks sur- 
rounding A stars and found that sandblasting has a neg- 
ligible effect on the observed structure, as radiation pres- 
sure blows out most incoming interstellar grains before 
they are allowed to intersect the disk. However, under 
this framework, only grains with /3 > 1 are ejected. Thus 
Figure [16] shows that radiation pressure does not protect 
the HD 61005 disk, as it does in A-stars. 

Nevertheless, even if radiation pressure does not pro- 
tect the disk against sandblasting, the stellar wind might, 
as only large interstellar grains with sizes greater than 
a few xO.l/im are allowed to enter astrosp heres freely 
without deflection (|Linde fc Gombosil 120001 ). Thus it is 
likely that interstellar sandblasting can only plausibly 
compete with interstellar gas drag if the astrosphere is 
smaller than or comparable to the observed debris disk 
size (Mann et al. 2006). The size of the HD 61005 astro- 
sphere is unconstrained by present observations. In gen- 
eral, observations and models of astrospheres surround- 
ing solar-type stars show sizes in the range ~ 10 — 10^ 
AU, depending on the ambient ISM and stellar wind con- 
ditions (see discussion in §5.1). Thus with a characteris- 
tic disk size of < 70 AU, it is not clear whether typical 
interstellar grains can intersect the HD 61005 disk. 

Detailed modeling of sandblasting is outside the scope 
of this paper. However, future theoretical work should 
investigate the effects of sandblasting on debris disks sur- 
rounding solar-type stars. Calculations of the ISM den- 
sity required for sandblasting to eject an observable flux 
of grains, the disk morphologies produced in this case, 
and the timescale for which sandblasting can be sus- 
tained would significantly aid in differentiating between 
this explanation and the gas drag models presented here. 

6. SUMMARY 

The morphology and polarization structure of HD 
61005 in the HST/ACS data (Figures [D - ^ strongly 
suggest that HD 61005 is a debris disk undergoing sig- 
nificant erosion by the ambient interstellar medium. The 
physical mechanism responsible for this erosion remains 
uncertain. Previous work has suggested that HD 61005 
may be interacting with an unusually dense cloud. How- 
ever, our high-resolution optical spectrum argues against 
this idea, instead suggesting an ambient ISM density typ- 



ical of local interstellar clouds. Thus the evolutionary 
state of HD 61005 may represent a commonplace, inter- 
mittent stage of debris disk evolution driven by interac- 
tion with typical, low-density gas. 

With this motivation, we considered the effects of secu- 
lar perturbations to grain orbits induced by ram pressure 
in warm, tenuous clouds. This mechanism can signifi- 
cantly distort grain orbits within a typical cloud cross- 
ing time and generate structures that very roughly repro- 
duce the HD 61005 images. Future work that incorpo- 
rates additional, more detailed physics may improve the 
agreement between the observations and interstellar gas 
drag models. The theoretical effects of interstellar sand- 
blasting for solar-type stars should also be investigated 
in greater detail. 

Regardless of the interpretation for HD 61005, we ex- 
pect interstellar gas drag is important at some level in 
shaping the structure and evolution of planetary debris 
disks. The frequency with which this effect is impor- 
tant strongly depends on the typical sizes, shapes, veloc- 
ities, and filling factors of warm interstellar clouds, which 
have poorly constrained global properties at present. 
Nevertheless, some morphological features common to 
nearby resolved debris disks (e.g., brightness asymme- 
tries, warps, and bow structures) can in principle be pro- 
duced in this way. A larger sample of spatially resolved 
debris disks at a wide range of wavelengths and more de- 
tailed theoretical work will help eliminate some of these 
remaining ambiguities. 
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APPENDIX 

RADIATION PRESSURE AND BLOW-OUT 

To provide a preliminary assessment of the unbound grain contribution, we compute the radiation pressure force 
according to lKruegell (|2003f )^Q: 

/2 
^{l-g^LoM^'F^du. (Al) 

Here, g^, is the grain scattering asymmetry, cjjy is the a lbedo, and Qff* is extinction cross section in units of the 
geometric cross section. H07 and ^Ca rpenter et al.l (|2008f ) showed that the stellar spectral energy distribution of HD 
61005 is well- matched by a main-sequence Teff = 5456 K Kurucz model atmosphere. We therefore use their best-fit 
spectrum in evaluating Eq. (|Aip to interpolate between the data and to extrapolate for the small fraction of missing 
flux longward of 24 /im and shortward of 0.3 /im. 

Figure [16] shows the ratio of the radiation pressure to gravitational force, P = Frp/Fq, for spherical particles 
with scattering properties computed using Mie theory. Results are shown for water ice (p = 1 g cm~^; Warren 

This expression is equivalent to that given bv lKohler Maiml (|2004l ) : 
Frp = / ^[Qf ^ + (1 - g.)Ql''']F^ du. 
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Fig. 16. — Ratio of the radiation pressure force to gravity (f3) for astrophysical silicate grains (left) and water ice (right) for the HD 61005 
system. The horizontal dotted line indicates the ratio above which grains launched by parent bodies on circular orbits become unbound. 
Thus silicate and ice grains in a conventional debris disk surrounding a solar type star are likely to remain bound to the star. On the other 
hand, the weak radiation field {(3 < 1) implies that radiation pressure alone does not impede interstellar grains from entering the system. 
Thus if the astrosphere surrounding HD 61005 is smaller than the disk, interstellar sandblasting could potentially erode the HD 61005 disk 
grains. 

1984) and Draine's astrophysical silicate (p = 3.5 g cm~^; Draine & Lee 1984). Mie theory radiation pressure 
calculations for silicate grains have been verified to within factors of a few using results from microwave analog 
laboratory data, computational Discrete Dipole Approximation (DDA) and T-matrix calculations, and solar system 
collection experiments (Wehry et al. 2004, Landgraf et al. 1999, and references therein). Solar system calculations 
performed on non-spherical, ballistic particle-cluster aggregates and ballistic cluster-cluster aggregates also yield similar 
results to spherical-grain Mie theory calculations (e.g.. Fig. 7 Mann et al. 2006, and references therein). For the various 
porosities shown in Figure [161 we used the Maxwell- Garnett rule to compute the approximate dielectric constant for 
a dilute medium. As applied to debr is disk systems, this ni ethod has been found to agree well with DDA calculations 
for aggregate porosities of P < 90% (jKohler Mannll2003 ). 

The calculation in Figure [TB] does not include the effect of stellar wind pressure, as the stellar wind parameters 
for HD 61005 are unknown, and the X-ray fiux of HD 61005 {Fx = 4 x 10^; Wichmann et al. 2003) exceeds the 
maxi mum value for which the Fx — relation of IWood et al.l (|2005l ) is reliable. Nevertheless, evaluating Equation 7 
from lStrubbe fc Chiang) (|2006l ) in the geometric optics limit, stellar wind pressure is predicted to be over an order of 
magnitude less significant than radiation pressure, even for mass loss rates approaching ~ 100 M©, approximately 
the highest mass loss rate yet observed. Thus the radiation pressure force of Figure [16] is expected to be representative 
of the total outward pressure force, independent of the stellar mass loss rate. 

If a grain is liberated from a parent body that is on a circular orbit, radiation pressure increases its eccentricity such 
that e = /3/(l — /?); for > 1/2, the grain is launched onto an unbound orbit. Debris disk grains tend to be highly 
porous with vacuum volume filling fractions of ~ 90% (e.g., Li & Greenberg 1998, Graham et al. 2007). Thus Figure 
[T6l shows that silicate and ice grains in a conventional debris disk surrounding a late-G star like HD 61005 are likely 
to remain bound to the star, unless the parent body population is highly eccentric. 

NEUTRAL GAS DRAG TOY MODELS 

Model Construction and Results 

The time for a 0.1 jiia grain to become unbound by interstellar gas drag is ~ 10^ — 10^ years, depending on the 
ambient ISM density, velocity, and relative flow direction, and the initial orbital elements of the grain (more precise 
estimates are given later in this section). This timescale may be compared to the collision time for submicron grains. 
Low eccentricity (e = 0.1 — 0.2) grains are pumped to moderate eccentricities (e = 0.3 — 0.4) b y interstellar gas drag 
within a typical orbital period of 500 yr. Thus, combining equations 7 and 9 from lChiang et all ([2009.) for moderately 
eccentric grains, the approximate collisional lifetime of submicron sized grains is 

2x10-^ H/R I^R/R f 0.95 MqY^^ f R Y^^fl-OA^^^ 

Here, the fractional luminosit y, L^^j /L^ = 2 x 10~ ^ , is ta ken from H07 and is identical to the result obtained from 
integrating the best fit SED in iRoccatagliata et al.l (|2009f ). Their fit implies that grains smaller than ^ 1 /xm are 
responsible for the bulk of the infrared excess and that these grains reside at a minimum distance of 96 ± 23 AU from 
the central star. This minimum disk radius is consistent with the characteristic radius of 70 AU adopted in Equation 

[HU 
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t = 500 yr t = 900 yr t = 1 600 yr t = 2800 yr t = 5000 yr 




Fig. 17. — Five snapshots of an initially unperturbed, near edge-on disk subjected to interstellar gas drag from a warm, low-density cloud. 
The model assumes a steady state develops such that each time a grain is lost from the initially unperturbed ring, a new replacement 
grain is generated within the initial "birth ring." The relative flow direction in this example is coplanar with the disk midplane, and the 
equilibrated disk is equivalent to the third column in the flrst row of Figure [12] The observed morphology does not change after 3 x 10^ 
years, setting the timescale for which the system is assumed to have achieved steady state. 

Since the collision time appears comparable to the time for the ISM to remove grains, we assume for our numerical 
models that a steady state is established in which the collisional generation of submicron grains within the posited 
birth ring is balanced by their removal by interstellar gas drag. We neglect the depletion of grains by collisions — 
this amounts to an order-unity error (see also Section 5.1, item 3). Within this framework, we follow the evolution 
of 10^ 0.1 jam particles. This single grain size is chosen for simplicity and is meant to be representative, given the 
polarization and color results described in §3. Based on the observed morphology, we choose an initial semimajor axis 
for all grains of 60 AU. Assuming the observed grains are collisionally created by grains on circular orbits. Appendix 
A suggests the initial grain eccentricities can range from e ~ to e ~ 0.4 for 0.1 /am grains, depending on the grain 
material and porosity. We adopt a representative value of e = 0.2 in our models, appropriate for moderately porous 
silicate grains. This eccentricity is consistent with the fractional ring width adopted in Equation IBll We further 
choose random inclinations drawn from a uniform distribution extending from 0° to 10°, based on an assumed disk 
aspect ratio of H/R ~ 0.1, consi stent with the system morphology in Figures [l| and [2] and with measurements of 
other highly inclined debris disks (jGolimowski et al]l2006l : [Kalasll2005l : iKrist et al.ll2005l ). Finally, we assume that the 
disk is initially circularly symmetric, such that the longitudes of ascending nodes, mean anomalies, and arguments of 
periastra are uniformly distributed between 0° and 360°. 

Having chosen initial orbital elements, we next calculate the secular perturbations to the angular momentum and 
Rung e-Lenz vectors from interste llar gas drag by averaging the Gaussian perturbation equations numerically over one 
orbit (jBrouwer fc Clemencelll96lh . We adopt an ISM density and encounter speed typical of warm interstellar clouds 
(nni = 0.2 cm~^, v^ei = 25 km s~^; Redfield 2006), though we note that the gas density is likely reduced by a factor a 
few inside the astrosphere relative to the nominal value outside (Bzowski et al. 2009). The direction of the flow vector 
relative to the disk has a signiflcant effect on the resulting disk morphology, and we therefore test a range of flow 
directions. 

To generate scattered-light images of the perturbed disk at some time after the initial encounter with the cloud, we 
calculate the classical elements from the resulting Runge-Lenz and angular momentum vectors and spread each grain 
over 100 points in its orbit in proportion to the time spent at each location. We estimate the equilibration timescale 
by monitoring the system morphology at successive times after the initial encounter with the cloud. Figure [iTl shows 
an example of the disk morphology evolution for a highly inclined disk, 10° from edge-on. The relative flow direction 
in this example is coplanar with the disk midplane. The observed morphology does not change after 3 x 10^ years, 
setting the timescale for which the system is assumed to have achieved steady state. 

To create scattered light model images for comparison to HD 61005, we incline the equilibrated disk by 80° and 
—80°. We also generate model images inclined by 0° (face-on), 45°, 90° (edge-on), and —45° to illustrate the more 
general consequences of neutral gas drag. In our notation, positive inclinations indicate that the lower half of the disk 
in the image is closer to the observer than the upper half. The pixel size in our images is 2.6 AU, corresponding to 
the projected size of a NICMOS pixel (0.0759'') at the distance of HD 61005 (34.5 pc). The scattered light images 
also assume a Henyey-Greenstein phase function with a scattering asymmetry typical of debris disks (g = 0.3; Meyer 
et al. 2007). If the grains in HD 61005 are similar to those in AU Mic, as suggested by the color and polarization 
results (§3), a larger scattering asymmetry may be more appropriate (e.g., g ~ 0.7; Graham et al. 2007). However, 
the adopted scattering asymmetry has a relatively small effect on the qualitative morphology for highly inclined outer 
disks like HD 61005, as the majority of scattering angles present are near 90°. 

The grid of models described above are shown in Figures [121 HSl and [TH The image display scale is logarithmic, 
and each model is ^310 AU (= 9'' for HD 61005) to a side. The labels at top indicate the disk inclination, and the 
labels at left indicate the azimuthal direction of the radial component of the flow vector. The bottom labels indicate 
the magnitude of the flow vector parallel and perpendicular to the disk midplane. The adopted coordinate system is 
illustrated with respect to the face-on disks in the left column of each Figure. 

When a signiflcant fraction of the relative flow vector is parallel to the disk midplane. Figures [12] and [13] show 
that a variety of morphologies can be produced, including bow struct ures (i = 0°), brightness asymmetries {i = 
80°, 90°, -80°), and warps {i = 90°/6> = 0°). As shown analytically bv iSchererl (|2000f ) and discussed in §4.2.2, the 
neutral gas drag force acts to rotate the grain pericenters into a direction perpendicular to the flow vector, resulting 
in a build up of particles in that direction. This effect is clearly observed in the face-on cases {i = 0°) in Figures [12] 
and [13] 

For a highly inclined disk and a small range of radial flow directions (0 ~ 0°, 180°), the models in Figures [12] and 
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[13] show a swept morphology somewhat similar to HD 61005. This effect is produced both by the build up of particles 
perpendicular to the flow vector and the non-zero inclination of the disk. Counterintuitively, the disk models with 
Vr = ^ and Vz = |'^| (Figure [14]) do not show this structure and instead produce a largely symmetric distribution about 
the disk midplane. As discussed in Scherer (2000), this is because the direction in which a given grain's pericenter 
rotates to become perpendicular to the flow depends on the initial pericenter direction. As a result, grains with an 
initial direction to pericenter above the nominal midplane will persist in having a direction to pericenter at or above 
the midplane as long as they are bound to the star. 
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